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1.  Introduction 


This  report  examines  basic  control  authority  characteristics  of  an  example  81 -mm  mortar 
equipped  with  a  controllable  internal  translating  mass.  Translating  mass  control  has  a  significant 
advantage  over  other  types  of  control  mechanisms,  such  as  canards  or  pulse  jets,  in  that  all 
components  of  the  control  mechanism  are  internal  to  the  projectile.  In  addition,  translating  mass 
actuators  can  be  implemented  using  simple  and  robust  hardware,  resulting  in  a  mechanism  less 
susceptible  to  damage  during  launch.  While  internal  translating  mass  control  has  been 
considered  for  many  types  of  projectiles,  its  use  on  mortars  is  particularly  attractive  due  to  a 
mortar’s  low  spin  rate  and  relatively  low  static  stability,  factors  that  reduce  actuator  power 
required  and  increase  control  authority. 

Control  authority  of  the  guided  mortar  is  examined  using  a  flight  simulation  of  an  example 
projectile  equipped  with  an  internal  translating  mass  (ITM).  This  simulation  comprises  a 
dynamic  model  of  the  projectile  coupled  to  a  control  system  for  the  ITM.  The  7-degree-of- 
freedom  dynamic  model  consists  of  14  scalar  differential  equations  that  are  integrated  forward  in 
time  from  a  given  set  of  initial  conditions.  The  ITM  control  system  utilizes  a  feedback 
linearization  controller  that  translates  the  mass  within  the  projectile  in  a  controlled  fashion  at  a 
frequency  equal  to  the  projectile  roll  rate.  While  an  outline  of  the  dynamic  model  is  presented 
here,  a  full  derivation  of  the  equations  of  motion  and  expressions  for  the  feedback  linearization 
controller  can  be  found  in  Rogers  and  Costello.1 

This  report  outlines  several  example  cases  and  trade  studies  to  demonstrate  the  capabilities  of 
translating  mass  control  with  the  81 -mm  mortar.  Configuration  parameters  such  as  mass  size, 
ITM  position  within  the  round,  and  ITM  oscillation  amplitude  are  all  varied  to  evaluate  their 
impact  on  system  performance.  Furthermore,  two  different  static  stability  configurations  are 
included  in  the  analysis  in  order  to  show  the  substantial  control  authority  improvements  when 
projectile  static  stability  is  reduced. 


2.  Dynamic  Model 


The  system  consists  of  two  major  components,  namely,  a  main  projectile  body  and  an  internal 
translating  mass.  The  main  projectile  body  is  largely  a  typical  mortar  projectile  with  the 
exception  of  an  internal  cavity  that  hosts  an  internal  mass.  The  internal  mass  is  free  to  translate 
within  the  main  projectile  cavity.  An  actuator  system  inside  the  projectile  exerts  a  force  on  the 


1  Rogers,  J.;  Costello,  M.  Control  Authority  of  a  Projectile  Equipped  With  an  Internal  Translating  Mass.  Journal  of 
Guidance,  Control,  and  Dynamics  2008,  31  (5),  1323-1333. 
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internal  mass  as  well  as  the  main  projectile  to  move  the  mass  inside  the  cavity  to  a  desired 
location.  A  schematic  of  the  81 -mm  mortar  projectile  equipped  with  an  ITM  is  shown  in 
figure  1,  and  figure  2  shows  the  earth-fixed  and  body-fixed  coordinate  systems. 


Figure  1.  Top  down  view  of  81-mm  mortar  equipped  with  an  ITM. 


Figure  2.  Position  coordinates  schematic  of  the  81-mm  mortar  projectile. 

Section  2.1  describes  the  dynamic  equations  of  the  ITM-mortar  system,  while  section  2.2 
outlines  the  physical  parameters  of  the  example  projectile.  Section  2.3  describes  the  ITM  control 
setup  and  configuration  parameters  used  throughout  the  results  section. 
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2.1  Dynamic  Equations 


The  dynamic  equations  for  this  7-degree-of-freedom  system  consist  of  a  set  of  kinematic, 
translational,  and  rotational  dynamic  equations.  The  reader  is  referred  to  Rogers  and  Costello1 
for  a  complete  derivation  of  the  dynamic  equations  for  a  projectile  equipped  with  an  internal 
translating  mass.  Only  the  final  results  of  the  derivation  are  summarized  in  this  report. 


Six  kinematic  equations,  relating  position  states  to  velocity  states,  are  identical  to  those  for  a 
6-degree-of-freedom  rigid  projectile  and  are  given  by  equations  1  and  2. 
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An  additional  kinematic  equation  is  given  by  the  trivial  relationship 

s  =  vs,  (3) 


which  relates  the  time  derivative  of  the  ITM  position  to  its  velocity  with  respect  to  the  projectile. 
Three  translational  dynamic  equations  are  also  identical  to  those  used  to  describe  rigid  projectile 
motion,  given  by 
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Another  translational  dynamic  equation  incorporates  the  acceleration  of  the  translating  mass  and 
is  given  by 
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Expressed  in  the  projectile  reference  frame,  the  components  of  the  rotational  kinetic  differential 
equations  are 
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The  dynamic  equations  of  motion  for  the  internal  translating  mass  projectile  are  collectively 
given  by  equations  1-5  and  10.  With  a  known  set  of  initial  conditions,  these  14  scalar  equations 
are  numerically  integrated  forward  in  time  using  a  fourth-order  Runge-Kutta  algorithm  to  obtain 
a  single  trajectory.  The  standard  aerodynamic  expansion  employed  for  projectile  flight  dynamic 
simulation  is  used  here,  with  aerodynamic  coefficients  generated  from  PROD  AS.  A  full 
description  of  the  weight  force  and  body  aerodynamic  forces  and  moments  is  provided  in  Rogers 
and  Costello.1 


2.2  Example  Projectile 

The  example  projectile  is  based  on  the  M821  81-mm  mortar.  Table  1  summarizes  dimensional 
and  inertial  properties  of  the  projectile. 

Table  1.  Properties  of  the  example  mortar  projectile. 


Diameter 

80.7  mm 

Projectile  Mass 

4.583  kg 

Axial  Inertia 

0.00371  kg-m2 

Transverse  Inertia 

0.04993  kg-m2 

CG  Stationline  (Nominal) 

0.282  m 

Fin  Cant  Angle 

1.4° 

Note:  CG  =  center  of  gravity. 


Note  that  the  CG  stationline  position  is  listed  as  “nominal,”  since  other  CG  positions  are 
considered  in  the  results  section  as  well.  Initial  conditions,  used  for  all  simulations  throughout 
the  results  section,  are  given  in  table  2. 


Table  2.  initial  conditions  for  the  example  mortar  projectile. 
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2.3  Internal  Translating  Mass  Control  Setup 

For  all  cases  considered  here,  the  cavity  that  contains  the  ITM  is  aligned  perpendicular  to  the 
projectile  axis  of  symmetry,  and  therefore  0T  =  0  and  y/T  =  nil  rad.  The  ITM  is  oscillated  at 

the  projectile  roll  frequency  according  to 


5  =  Acos(p)  ,  (14) 

where  s  is  the  ITM  distance  from  the  center  of  the  cavity,  A  is  the  ITM  oscillation  amplitude 
(exactly  half  of  the  ITM  stroke),  and  (p  is  the  projectile  roll  angle.  Figure  3  shows  a  diagram  of 
how  translating  mass  motion  occurs  through  different  portions  of  the  projectile  roll  cycle. 


Figure  3.  Moving  mass  position,  velocity  ( V  )  and  acceleration  ( a  )  over  one  roll  cycle,  viewed  from 
the  rear  of  the  projectile.  Mass  motion  is  sinusoidal  and  has  the  same  frequency  as  the 
projectile  roll  rate. 


ITM  oscillations  are  initiated  as  soon  as  the  projectile  is  launched  and  result  in  full  control 
authority  in  the  inertial  frame  negative  y  direction  (to  the  left  when  looking  downrange).  Two 
sets  of  ITM  configuration  parameters,  given  in  table  3,  are  studied  in  order  to  detennine  how 
different  variables  affect  control  authority. 


Table  3.  ITM  configuration  parameters  for  example  trajectories. 


Projectile  Physical  Variable 

Parameter  Set  A 

Parameter  Set  B 

ITM  mass  size  (g) 

125 

125 

ITM  stroke  (cm) 

1.27 

5.08 

Projectile  CG  stationline  (m  from  base) 

0.282 

0.209 

Cavity  offset  from  CG  (m) 

0.167 

0.167 

Parameter  set  A  represents  the  configuration  originally  proposed  by  the  U.S.  Army  Research 
Laboratory,  while  parameter  set  B  represents  a  projectile  with  reduced  stability  (CG  farther  aft) 
and  a  significant  increase  in  ITM  oscillation  amplitude  (“stroke”).  A  stroke  of  5.08  cm  was 
determined  to  be  approximately  the  maximum  amplitude  possible  given  the  width  of  the  round  at 
that  cavity  offset.  Also,  note  that  a  CG  stationline  of  0.209  m,  used  in  parameter  set  B,  ensures  a 
static  margin  of  at  least  0.25  cal.  throughout  the  flight. 
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3.  Results 


Both  example  cases  and  trade  studies  were  perfonned  to  demonstrate  the  capabilities  of  ITM 
control.  Section  3.1  outlines  example  controlled  simulations  using  both  sets  of  configuration 
parameters,  as  well  as  uncontrolled  cases  in  order  to  establish  a  baseline  for  comparison.  Section 
3.2  includes  trade  studies  that  examine  the  effect  of  ITM  size,  stroke,  and  cavity  offset  on  control 
authority. 

3,1  Example  Trajectories 

Several  controlled  trajectories  are  simulated  and  compared  to  uncontrolled  simulations  (in  which 
the  mass  is  present  but  not  oscillated).  As  described  in  section  2.3,  the  controller  commands 
mass  oscillation  to  begin  immediately  after  launch  and  continues  oscillation  throughout  the 
remainder  of  flight.  The  oscillation  function  given  in  equation  14  corresponds  to  control  applied 
in  the  inertial  frame  negative  v  direction  (see  figure  2).  Furthermore,  the  set  of  initial  conditions 
given  in  table  2  corresponds  to  approximately  maximum  range  for  the  example  mortar  projectile. 
Therefore,  all  controlled  trajectories  show  maximum  control  authority  in  deflection  for  a  specific 
ITM  configuration  at  maximum  range.  This  maximum  deflection  capability  is  used  as  the 
quantitative  measure  of  control  authority. 

Figures  4-10  show  four  example  trajectories.  Two  uncontrolled  trajectories  are  shown,  one  for  a 
projectile  of  nominal  stability  (CG  position  of  0.282  m)  and  one  for  a  projectile  of  reduced 
stability  (CG  position  of  0.209  m).  Note  that  these  values  correspond  to  those  used  in  ITM 
parameter  sets  A  and  B,  respectively.  Several  states  of  the  reduced  stability  uncontrolled 
trajectory  have  been  excluded  since  they  are  nearly  identical  to  those  of  the  nominal  stability 
uncontrolled  case.  Two  controlled  trajectories  are  also  shown,  one  using  parameter  set  A  and  the 
other  using  parameter  set  B. 

The  most  significant  plots  are  figures  5,  deflection  vs.  range,  and  10,  angle  of  attack  vs.  time. 
Both  show  that  the  ITM  control  mechanism  has  little  effect  for  the  configuration  given  by 
parameter  set  A.  Maximum  control  authority  for  this  case  yields  deflection  of  ~1  m  from 
nominal  at  maximum  range  and  a  trim  angle  of  attack  of  <0. 1°.  However,  parameter  set  B  shows 
significantly  better  results,  with  a  deflection  of  over  20  m  from  nominal  at  maximum  range  and  a 
trim  angle  of  attack  of  ~0.3°.  Figure  5  also  shows  trajectory  deflections  for  the  uncontrolled 
rounds.  This  is  due  to  drift,  which  is  caused  by  an  interaction  between  the  projectile  rolling  and 
pitching  dynamics  with  the  gravity-induced  trajectory  curvature.  Drift  is  larger  for  the  reduced 
stability  configuration  than  for  the  nominal  stability  case  due  to  a  decrease  in  restoring 
aerodynamic  moment. 

Figures  5  and  10  highlight  the  fact  that  ITM  systems  typically  only  show  good  perfonnance  in 
rounds  with  reduced  static  stability,  where  ITM  size  and  stroke  are  as  large  as  possible. 
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Range  (m) 


Figure  4.  Altitude  vs.  range.  All  trajectories  overlay  one  another. 


Range  (m) 


Figure  5.  Deflection  vs.  range. 
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Figure  6.  Total  velocity  vs.  time.  All  trajectories  overlay  one  another. 


Figure  7.  Roll  rate  vs.  time.  All  trajectories  approximately  overlay  one  another. 
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Time  (s) 


Figure  8.  Zoom  view  of  roll  rate  vs.  time.  Notice  roll  rate  oscillations  as  ITM  oscillates, 
changing  the  axial  moment  of  inertia  of  the  system. 
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Figure  9.  Segment  of  ITM  displacement  from  centerline  vs.  time. 
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Figure  10.  Angle  of  attack  vs.  time. 

3.2  Trade  Studies 

A  trade  study  was  conducted  to  examine  the  effect  of  cavity  offset  on  ITM  control  authority.  In 
addition,  the  ITM  stroke  was  varied  for  each  cavity  offset.  Simulations  utilizing  identical  initial 
conditions  to  the  examples  in  section  3.1  were  performed  using  both  nominal  stability  (CG 
stationline  of  0.282  m)  and  reduced  stability  (CG  stationline  of  0.209  m).  Radial  deviation  from 
the  nominal  case  was  determined  for  each  simulation.  Note  that  radial  deviation  is  calculated  as 
the  square  root  of  the  sum  of  the  squares  of  the  deviations  in  deflection  and  range,  which  occurs 
due  to  a  slight  off-axis  response  of  the  projectile.  Figure  1 1  shows  control  authority  results  for 
the  nominal  stability  cases,  while  figure  12  shows  control  authority  results  for  the  reduced 
stability  cases.  ITM  size  was  125  g  for  all  simulations.  Note  that  deviation  is  nearly 
independent  of  cavity  offset.  This  is  due  to  the  fact  that  the  control  mechanism  relies  on  a  lateral 
mass  center  offset,  which  can  be  produced  at  any  stationline  on  the  projectile.  Slight  increases  in 
radial  deviation  seen  with  larger  offsets  are  the  result  of  drag  penalties  associated  with  higher 
angle  of  attack  oscillations  at  large  cavity  offsets.  Also  note  that  not  all  combinations  of  ITM 
stroke  and  cavity  offset  are  physically  feasible  given  the  constraints  of  the  81 -mm  round. 
However,  the  results  show  important  trends,  specifically  that  only  large  ITM  stroke  combined 
with  reduced  stability  can  produce  any  useful  control  authority. 
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Figure  11.  Radial  deviation  from  standard  trajectory  vs.  cavity  offset  for  nominal  stability  (ITM 
mass  of  125  g,  CG  stationline  of  0.282  m). 
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Figure  12.  Radial  deviation  from  standard  trajectory  vs.  cavity  offset  for  reduced  stability  (ITM  mass 
of  125  g,  CG  stationline  of  0.209  m). 
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Similar  trade  studies  were  conducted  by  fixing  the  cavity  offset  at  0. 167  m  and  varying  the  ITM 
stroke  and  mass  size.  Again,  the  study  was  conducted  for  both  the  nominal  and  reduced  stability 
cases.  Figure  13  shows  control  authority  results  for  the  nominal  stability  cases,  while  figure  14 
shows  control  authority  results  for  the  reduced  stability  cases. 


Figure  13.  Radial  deviation  from  standard  trajectory  vs.  ITM  stroke  for  nominal  stability  (cavity 
offset  of  0. 167  m,  CG  stationline  of  0.282  m). 

Figures  13  and  14  again  demonstrate  that  the  reduced  stability  ITM  projectile  exhibits 
significantly  better  performance.  A  larger  stroke  and  larger  mass  produce  much  greater  control 
authority.  Flowever,  it  is  important  to  note  that  larger  mass  sizes  and  oscillation  amplitudes 
require  larger  actuators  and  more  battery  power. 


4.  Conclusion 


A  control  authority  study  of  an  81 -mm  mortar  equipped  with  a  controllable  internal  translating 
mass  was  presented.  Example  results  show  that  reasonable  control  authority  is  achieved  only  in 
the  case  of  decreased  static  stability.  Parametric  trade  studies  further  showed  that  control 
authority  increases  as  mass  size  and  oscillation  amplitude  increase,  and  is  insensitive  to  cavity 
offset. 
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Figure  14.  Radial  deviation  from  standard  trajectory  vs.  ITM  stroke  for  reduced  stability  (cavity 
offset  of  0.167  m,  CG  stationline  of  0.209  m). 

Analysis  of  these  results  leads  to  the  conclusion  that  use  of  an  ITM  mechanism  results  in  far  less 
control  authority  than  is  possible  with  a  mechanism  such  as  canards,  which  can  achieve  hundreds 
of  meters  of  swerve  in  deflection  at  maximum  range.  Furthermore,  in  order  to  create  reasonable 
control  authority  with  the  ITM  mechanism,  the  projectile  static  margin  would  need  to  be 
substantially  reduced  to  -0.25  cal.  or  less.  Therefore,  since  the  goal  of  controlled  munitions  is 
often  primarily  to  reduce  or  eliminate  trajectory  dispersion,  it  is  highly  unlikely  that  an  81 -mm 
mortar  equipped  with  an  ITM  control  system  would  be  capable  of  control  maneuvers  needed  to 
eliminate  launch  errors  and  trajectory  disturbances. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


cy  The  cosine  of  argument  angle  ? . 

CG  Center  of  gravity. 

f input  Scalar  value  of  the  input  force  exerted  by  the  controller. 

g  The  acceleration  due  to  gravity  (32.2  ft/s"). 

IP  Mass  moment  of  inertia  matrix  of  the  projectile  body  with  respect  to  the  projectile 

reference  frame. 

IT  Mass  moment  of  inertia  matrix  of  the  internal  translating  mass  with  respect  to  the 

projectile  reference  frame. 

ITM  Internal  translating  mass. 

1 1,  J  i,  K  i  Inertial  reference  frame  unit  vectors. 

I P,J P,K p  Projectile  reference  frame  unit  vectors. 

IT  ,JT  ,Kt  Internal  translating  mass  reference  frame  unit  vectors. 


mP  Projectile  body  mass. 

mT  Internal  translating  mass. 

m  Total  system  mass. 

M  x ,  MY ,  My  External  moment  components  on  the  proj  ectile  body  expressed  in  the  proj  ectile 
reference  frame. 

PMP'  Components  of  the  angular  velocity  vector  of  the  projectile  body  expressed  in  the 

projectile  reference  frame. 

nc^T  Distance  vector  from  the  center  of  mass  of  the  system  to  the  internal  translating 

mass  center  of  mass. 
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rp^T  Distance  vector  from  the  projectile  center  of  mass  to  the  internal  translating  mass 

center  of  mass. 

s  Position  of  the  internal  translating  mass  along  its  line  of  movement  with  respect  to 

the  center  of  the  cavity. 

sr  The  sine  of  argument  angle  y . 


Transformation  matrix  from  the  projectile  reference  frame  to  the  translating  mass 
frame. 


t  The  tangent  of  argument  angle  y . 

u,v,w  Translation  velocity  components  of  the  composite  body  center  of  mass  resolved  in 

the  projectile  reference  frame. 

vc/1  Velocity  of  the  system  mass  center  with  respect  to  the  inertial  frame. 


vs  Magnitude  of  the  velocity  of  the  translating  mass  with  respect  to  the  translating 

mass  reference  frame. 

X’T’Z  Position  vector  components  of  the  composite  body  center  of  mass  expressed  in  the 

inertial  reference  frame. 

X,Y,Z  Total  external  force  components  on  the  composite  body  expressed  in  the 
projectile  reference  frame. 

9,  V  Euler  roll,  pitch,  and  yaw  angles. 

0T,  y/T  Euler  pitch  and  yaw  angles  for  the  orientation  of  the  line  of 

movement  of  the  internal  translating  mass  with  respect  to  the 
projectile  body. 

(bpn  Angular  velocity  of  the  projectile  body  with  respect  to  the  inertial  frame. 

CA,  ( j)  Vector  component  operator  that  outputs  a  column  vector  comprising  the 

components  of  the  input  vector  y  expressed  in  reference  frame  X . 

( j)  Cross  product  operator  that  outputs  a  skew-symmetric  matrix  using  the 

components  of  the  input  vector  y  expressed  in  reference  frame  X . 
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